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Abstract 

We discuss an interesting possibility of nontrivial, quark-mass induced chiral 

vacuum alignment in color-flavor locking phase of cold, dense QCD. With 

the simplifying assumption that the gaps for quarks are identical to those of 

antiquarks and the light quark masses are given by m u = nid and m s /md = 15, 

we find the true chiral vacuum can align only to one of the discrete number 

of directions in the continuum of chiral vacua. The alignment depends on 

the size of the diquark condensates, and the vacuum transitions between the 

discrete vacua caused by the evolution of the diquark condensates can be first 

order phase transition with vanishing or nonvanishing latent heat, depending 

on the vacua involved. It is also shown that as /i — > oo, where fi is the 

baryon chemical potential, parity is spontaneously broken through the vacuum 

alignment. 
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The chiral symmetry SUl(3) x SUr(3) of quantum chromodynamics (QCD) with three 
massless quark flavors is spontaneously broken to SUl + #(3) at low energies. Any point in the 
continuum of the chiral vacua, the coset space [SUl(3) x SUr(3)]/SUl + .r(3), can be chosen 
as the vacuum, since any two vacua are equivalent as far as physics is concerned. However 
when quarks receive small (current) masses the chiral symmetry is then explicitly broken, 
and the degeneracy of the vacua is lifted. The quark mass term then picks up the lowest 
energy state in the continuum of the vacua as the true chiral vacuum. This is called Dashen's 
chiral vacuum alignment To determine the true vacuum one has to find the quark- mass 
induced potential that lifts the degeneracy of the vacua. Since the Nambu-Goldstone bosons 
of the chiral symmetry breaking are small fluctuations about a point in the continuum of 
the vacua, the potential is nothing but the meson mass term in the chiral Lagrangian of the 
nonlinearly realized Nambu-Goldstone bosons. Therefore, the true vacuum is the one that 
minimizes the potential 

K(S) = -£ m (S) oc -Re[Tr(mS)], (1) 

where S G SU(3) denotes the chiral fields for the Nambu-Goldstone bosons, while £ m (S) is 
the quark-mass induced meson mass term in the chiral Lagrangian. 

When the true vacuum is nontrivial, that is, So 7^ I, where So is the vacuum that 
minimizes the potential, some interesting phenomena could arise. For instance, if the quark 
mass matrix were given as m oc Diag(— 1, —1, —5) with S > 1/2, then the vacuum S would 
have an imaginary part in its matrix elements 0. Then writing S = S exp(iH A T A ) , where 
U A and T A denote the octet pseudo scalar mesons and the SU(3) generators, respectively, 
it can be easily seen that the mass term ([I]) could give rise to a spontaneous CP violation, 
which enables, for example, rj decay into two pions. Of course, with the quark mass matrix 
realized in nature, the true vacuum is at S = /, and there is no such CP violation. 

Recently, cold, dense quark system with large baryon chemical potential fi has received 
strong interest |2|||. As well known, the Fermi surface of such a dense system is unstable 
against Cooper pairing when an attractive force is present. At large /i a dressed gluon 
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exchange (in hard dense loop approximation) provides such a force, and causes the system 
to be in color-flavor locking phase in which quarks condensate in a pattern 

ex -[fe + A; 2 (Vb) ,-^) 6 ,], (2) 

where Xi\ V 9 ^ with a the color index, z = 1, • • • , 3, the flavor index, denote the two-component 
Weyl fermions for the left-handed quarks and the complex conjugate of the right-handed 
quarks, respectively. The unitary matrices Uq, Vq £ U(3) can be arbitrary in the absence of 
axial U(l) anomaly, but the anomaly effect, even though small because of the suppression 
of instanton effects at high density [f§,|5|, chooses a parity even vacuum ]5]|| in which 

U = V . (3) 

Upon the condensation of the quarks the symmetry of dense, massless QCD, SUl(3) x 
SUr(3) x Ux(l) x XJ B (1), where Ua(1) and Ub(1) denote the axial and the baryon num- 
ber symmetry, respectively, is spontaneously broken to SUi + ij(3), generating 10 Nambu- 
Goldstone bosons (mesons). The Ua(1) is not an exact symmetry, but at large /i, due to 
the suppression of anomaly, can be regarded a good approximate symmetry. As in vacuum 
QCD a continuum of chiral vacua [SUr,(3) x SU#(3) x Ua(1) x Ub(1)]/SUl+,r(3) arises upon 
the spontaneous symmetry breaking. Note, however, that the vacua connected by the Ua(1) 
rotation are only approximately degenerate, because the Ua(1) is not an exact symmetry. 
It turns out that the vacuum of lowest energy is parity even, hence the relation (|3]). 

Now, when the current quark masses are turned on, the degeneracy of the chiral vacua 
is lifted. The true vacuum can be picked up, as in vacuum QCD, by minimizing the quark 
mass induced potential V(E) = — £ m (E) of dense QCD. 

The meson mass term £ m (£) in color-flavor locking phase has a different form than 
in vacuum QCD because of the absence of left-right quark condensates. For small quark 
masses, £ m (E) can be expanded in powers of the quark mass matrix. The absence of a left- 
right quark condensate due to the suppression of instanton effects renders the leading term 
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to be quadratic in quark mass @. The most general form for the leading £ m (E), consistent 
with the chiral symmetry and the condensates (0), is given by |7|-p 

£ m (£) = A [Tr(m*£)f + B Tr[(m*S) 2 ] + C Tr^S)^^*^ 1 ) + H.c. , (4) 

where S G U(3) denotes the chiral fields for the 9 mesons. Note that the Nambu-Goldstone 
boson associated with the baryon number symmetry remains exactly massless, and thus does 
not appear in the mass term. 

The coefficients A, B and C can be determined either by matching the vacuum energy 
of (|4]) with that computed in the microscopic theory || or by integrating out quark fields 
in the effective Lagrangian of quarks and the mesons 0. Here we take the latter approach. 
Using the global color-chiral-axial-baryon number symmetry of dense QCD and the pattern 
of the diquark condensates (0) one can write an effective chiral Lagrangian for the quarks 
and the Nambu-Goldstone bosons as 

C = i%t° V d»xt + VXt*\ a i + m^dvtf - fM<p?d°<pt - [m ijX ?tf + H.c] 



+ 



X?(At )f jX ) + tfCA^} + H.c] + C NG (U, V), (5) 



where 



(At)g- = hU* a U* b + k 2 U* a U* b , 

{Av)$ = -[hiV?V} +hVfV*}. (6) 

The unitary matrices U and V denote the nonlinearly realized Nambu-Goldstone bosons 
arising from the symmetry breaking SU C (3) x SUl(3) X Ua+b(1) — > SU c+ l(3) through the 
XX condensation and SU C (3) x SUr(3) x Ub-a{^) ~^ SU c+ r(3) through the (p(p condensa- 
tion, respectively. The C NG (U,V) is the usual chiral Lagrangian for the Nambu-Goldstone 
bosons alone [0,0 • We also note that with (^) it was assumed, for simplicity, the gaps for 
antiquarks are identical to those for quarks . It can be easily seen that the most important 
interactions for our calculation, those between the quarks and the Nambu-Goldstone bosons, 
satisfy the Goldberger-Treiman relation. This Lagrangian can be regarded as an effective 



Lagrangian before color is gauged for the quarks and the 18 Nambu-Goldstone bosons U, V. 
Upon gauging color it can be seen that 8 out of the 18 Nambu-Goldstone bosons are eaten 
by the gluons via Higgs mechanism and there remain 10 color-singlet mesons as low energy 
excitations. 

Integrating out the quark fields, which corresponds to the evaluation of quark one-loop 
diagrams in the background of constant U, V fields with two Dirac mass insertions (see Fig. 
1), we obtain the mass term (H) in which E is given by 

E = UV ] (7) 

and 

A = ij {(*? + kl)h-{p)h + ip) + h(h + k 2 /3)[I 8 -(p)h + (p) 
-/i+(p)/ 8 _(p)] + h(h + fe 2 /3)[/ 8 _(p)/ 8+ (p) - /!_(p)/ 8+ (p)] 
+ k 2 /3) 2 [I s 4p)I 8+ (p) - h-{p)h+(p) - h + (p)h-(P) 

+/l_(p)J 1+ (p)]}, 

B = iJ -0^ {2W 8 -(p)/ 8+ (p) + k 2 {k x + fe 2 /3)[/ 8 _(p)J 8+ (p) 

-J 1+ (p)/ 8 _(p)] + fc 2 (*i + fc 2 /3)[J 8 _(p)/ 8+ (p) - J 1 _(p)/ 8+ (p)]} , 

-/ 1+ (p)J 8 _(p) + / 1 _(p)J 1+ (p)}, (8) 

where 

J lT (p) = 1/hpg + (bl T + m?], J 8T (p) = 1/hpg + (bl T + mg]. (9) 

Here mi and m§ are the Majorana masses for the singlet and octet quarks, respectively, 
under the unbroken SUl+j?(3) and are given as 

ml = (3&i + k 2 ) 2 , mj = kl (10) 
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FIG. 1. Diagrams to be evaluated in the constant background of the chiral fields U and V. 

Generally the diquark condensates depend on energy, but here we shall ignore this fact 
and treat ki as constants. Then the integration over the loop momentum can be easily 
performed by doing contour integration over p Q first and then replacing d 3 p — > 4irfi 2 d\p\ 
for the integration over the spatial components. It can be easily seen that A and B arise 
from the first diagram in Fig.l with A, B cx A 2 ln(/z 2 /A 2 ), while C, coming from the second 
diagram, is given by C oc A 4 /// 2 ln(/z 2 /A 2 ), where A ~ k{. Note that C is suppressed by a 
factor A 2 //z 2 compared to A, B. 

To determine the true chiral vacuum at a given chemical potential we have to minimize 
the potential V(E) = — £ m (S), with £ m (£) given by @ and (§). To demonstrate a non- 
trivial vacuum alignment, we shall now take a simplified form for the quark mass matrix, in 
which the up and down quark masses are identical, as 



With this quark mass matrix it is easy to see that the E that minimizes the potential 
must be of the form 



m = m d Diag(l, 1,5), 5 = m s /m d . 



(11) 




(12) 



where a, (3 are complex variables. 



Substituting (fLip and (|T^) into (|J) we obtain 



V = -2m 2 d {Re[A(2a + (3d) 2 + B(2a 2 + (3 2 6 2 )} + C \2a + (35\ 2 } , 



(13) 



which can also be written as 
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V = -8m 2 d {{2A + B) cos 2 6 + {A + B)5 2 /2 cos 2 + A5 cos(fl + 0) + C5 cos(# - 0) } (14) 

by putting a = exp(i9) and (3 = exp(i0). 

We first notice that the potential is symmetric under the transformation (at, (3) — > 
(—a, — so the minima of the potential must occur in pairs. Secondly we observe that for 
arbitrary coefficients A, B, and C the potential is stationary when 9, satisfy 

sin(# ± 0) = 0, sin(2fl) = 0, sin(20) = 0, (15) 

which have 8 common solutions corresponding to the following (a,/3) pairs 

(-i, -i), (-1, 1), (1, 1), (i, -i) (16) 

and their partners of opposite sign. Of course none of these pairs needs necessarily minimize 
the potential, but it can be shown numerically that the minima occur always on one of 
these solutions. This shows that the true chiral vacuum can align only to one of these eight 
directions. 

In ideal situation we may know the dependence of the gap parameters hi on /x, could 
choose the true chiral vacuum at a given chemical potential from the above discrete vacua, 
and investigate vacuum transitions as the chemical potential evolves. However, presently 
there is no reliable calculation to determine the gaps as functions of the chemical potential 



except when the chemical potential is extremely large fL4 |, in which case Shwinger-Dyson 



equation can be used as an approximate gap equation and solved fl5|-|17|. Even in this 
case the absolute magnitudes of the gaps are yet to be determined, but, it may not be so 
unreasonable to assume that hij \x are in the order of 0.1 at large chemical potential P,[16|,|6| . 

Taking into account this uncertainty we shall here treat hi as free parameters and study 
the vacuum transitions as hi vary. For definiteness, we shall put 5 = 15 and scan numerically 
the minima of the potential in the hi parameter space defined by < k\f 'fi < 0.5 and 
—0.5 < hi) 'fi < 0.5. Note that hi can always be assumed positive, since its phase can be 
rotated away by the baryon number symmetry of the Lagrangian (|5|). 
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The result of the numerical scanning is shown in Fig. 2. As we see, the parameter 
space is divided into four domains according to their vacuum directions. In this figure the 
corresponding vacua of opposite sign to those in the group ( |16D are not included. The reason 
for this is that the numerical scanning shows there is always a potential barrier between a 
vacuum in the group ( Pp and the partners of those in the group ( p.6[ ) , so in infinite volume 
limit the transitions between the vacua fll6|) and their partners are negligibly small and 
can be ignored. Therefore, only the transitions within the group of vacua in ([!(]) need be 
considered. 




0.1 0.2 0.3 0.4 



k 1 

FIG. 2. Gap parameter space is divided into domains according to their vacua. Each domain 
is associated with a unique vacuum. k{ are in the unit of fj, = 1. 

Because the vacuum can align only to a discrete number of directions the transitions 
between vacua will be of first order phase transition. To understand what happens at the 
transition, we look more carefully at the potential on the boundaries between the domains. 
For convenience, we shall designate the domains associated with the vacua defined in ( |T6|) 
by (I), (II), (III), and (IV), respectively. 

First, consider the transition between the domain (I) and (II). As one approaches the 
boundary from either side of (I) and (II) it can be shown numerically that a potential valley 
opens up in (6, 0) space along the line = —9 — tt, with the bottom of the valley connecting 
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the two points (— tt/2, — tt/2) and (— 7r,0). When exactly on the boundary, we have from 
V(£ r ) = V(E n ), with S MI defined by (0) and ([16]), 

(2A + 5) + (A + B)5 2 /2 - 2C5 = 0. (17) 

Then it is easy to see that on the bottom of the valley (i.e. along the direction <fi = —9 — tt) 
the potential is constant with V = &m 2 d (A — C)S. Thus, in this case there will be no latent 
heat released at the phase transition. 

Similarly, for the transition between (II) and (III) it can be shown that a valley opens 
up along the line = as the boundary is approached, and that the relation 

A + C = 0, (18) 

holds on the boundary. On the bottom of the valley the potential is given by 

V = -8m 2 d {{2A + B) cos 2 6 + {A + B)5 2 /2} , (19) 

which shows a barrier between the two vacua. This then indicates that there will be latent 
heat released at the phase transition. This may have some important phenomenological 
consequence in dense systems. 

In a similar fashion we can easily show that for the transition across the domain (III) 
and (IV) a potential valley opens up along the direction = —9, and on the bottom of the 
valley the potential is constant with V = —%m 2 d (A — C)5. In this case there will be no latent 
heat released as in the transition between (I) and (II). Also for the transition between (I) 
and (IV) a valley opens up along = —tt/2 direction and on the bottom of the valley the 
potential is given by 

V = -4m 2 {A + B)5 2 cos 2 e, (20) 

which shows a barrier between the two vacua, and consequently, nonzero latent heat at the 
transition. 

Could the phase transition of this kind occur when the chemical potential increases from 
zero to an asymptotic value? Although it is difficult to answer this question conclusively until 
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we have the quark-mass induced meson potential at an arbitrary value of /i, which at smaller 
H would probably contain the instanton-induced potential that gives mass to the rj' meson 
and the old Tr(mS) as well as the quark- mass quadratic £ m (S) in (|]), there is an interesting 
observation concerning this question. It is well known that at large chemical potential the 
sextet components of the condensates (fj) is suppressed |T%|JEj|] , thus ki ~ — k 2 , which then 
suggests the system must be in domain (I) at high density. The vacuum associated with 
domain (I) is So = Diag(— i, —i, —1) that has an overall factor —% compared to the vacuum 
Sq = I at zero density. Although this does not imply a first order phase transition of the 
kind hitherto discussed, at least it does suggest that the vacuum must shift from the unit 
matrix at zero density to something else as the chemical potential increases. 

This also leads to an interesting consequence of the vacuum alignment, namely that 
parity must be spontaneously broken at high density. The possibility of parity violation in 
dense QCD was pointed out by Pisarski and Rischke 0], and noted also in pO|j21| P1, who 



observed that the axial U(l) is almost a good symmetry at large \x due to the suppression of 
instanton effects, so a vacuum with no parity symmetry can be as easily formed as the parity 
even vacuum. Since the axial U(l) is never exactly restored this possibility, however, cannot 
be realized with massless quarks. The lowest energy state is always the parity even vacuum. 
However, when quarks are massive, the explicit breaking of the axial U(l) by quark masses 
can be more important than the small anomalous breaking, so the true vacuum could break 
parity. 

Indeed, we see that this possibility is realized through the chiral vacuum alignment. Since 
the diquark condensates satisfy (|3]) in parity even vacuum it can be seen using (0) that only 
the vacuum So = / associated with the domain (III) is parity even. Therefore, when the 
system is in domain (I), (II), and (IV) parity is spontaneously broken. Moreover, since the 
system should be in domain (I) as fi — > 00, parity must be broken at high density. 

To conclude, we have studied quark-mass induced chiral vacuum alignment in cold, dense 
QCD. When quarks are massless any point in the continuum of chiral vacua can be chosen 
as the vacuum, but when nonzero quark masses are introduced, the true vacuum must be 
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found via Dashen's procedure. We have shown that in color-flavor locking phase at large 
chemical potential the vacuum can align only to a discrete number of directions, and the 
nature of the vacuum transitions is of first order with vanishing or nonvanishing latent heat. 
It was also shown that at high density parity is spontaneously broken through the vacuum 
alignment. The consequence of the first order phase transitions and the parity violation may 
become important in dense systems like quark stars and heavy ion collisions. 

Finally, we remark on the calculation of the meson masses. Usually in meson mass 
calculation it was assumed that the vacuum is at S = /, and the meson potential was 
expanded around this vacuum to pick up the meson spectrum. However, as we have seen, 
the vacuum is not necessarily at the unit matrix and so this is not always correct. For 
correct meson masses one should first find the true vacuum and then expand the potential 
about it. 

Useful comments by D.K. Hong are acknowledged. This work was supported in part by 
the Korea Science and Engineering Foundation (KOSEF). 
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